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Abstract 

We have searched with the KLOE detector for the P and CP violating decay 77 — > 
it + it~ in a sample of 1.55 x 10 7 r/'s from the decay <\> — > 777 of ^-mesons produced 
in e + e~ annihilations at DA<1?NE. No signal is found. We obtain the upper limit 
BR(r? -> 7r + 7r~)< 1.3xl0~ 5 at 90% confidence level. 
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The study of 77 decays provides an excellent laboratory for testing the validity 
of symmetries of the physical world. The decay 77 — > 7r + 7r~ violates both P and 
CP invariance. In the Standard Model, the reaction can proceed only via the 
weak interaction with a branching ratio of order 10 -27 according to Ref. [1]. 
Higher branching ratios are conceivable either introducing a CP violation term 
in the QCD lagrangian through the so-called 9 term (a branching ratio up to 
10~ 17 can be obtained in this scheme compatible with the experimental limit 
on the neutron electric dipole momentum) or allowing a CP violation in the 
extended Higgs sector (in this case 10 -15 can be reached) as described in 
Ref. [1]. The detection of this decay at any level accessible today would signal 
P and CP violation from new sources, beyond any considerable extension of 
the Standard Model. The best published previous result from a direct search 
for the 77 decay to 7r + 7r~ was obtained in 1999: BR(r) — > 7r + 7r~) < 3.3 x 10~ 4 
at 90% confidence level [2]. 

We present here the result of a direct search for this decay performed with the 
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KLOE experiment at the Frascati factory, DA$NE, based on an integrated 
luminosity of 350 pb" 1 collected in the years 2001 and 2002. DA$NE is an e + e~ 
collider working at the centre of mass energy of 1019.5 MeV, i.e., at the mass. 
mesons are produced nearly at rest in the laboratory. Beams in DA$NE 
collide with a crossing angle of ir — 0.025 rad. 0-mesons are therefore produced 
with a momentum \p^,\ ~ 12.5 MeV/c in the horizontal plane, directed towards 
the centre of the storage rings. The precise value of the momentum, together 
with the centre of mass energy and the beam spot position, is determined run 
by run using Bhabha scattering events. 

The total cross section for e + e~^0 is ~3 fib. r\ mesons are copiously produced 
at a typical rate of ~2 Hz through the radiative decay — > 777, which has a 
branching ratio of 1.3%. In the decay chain — > 777 — > 7r + 7r~7, the photon 
is monochromatic (i? 7 = 363 MeV) and is emitted according to an angular 
distribution dN '/ 'dcosO^ ~ (1 + cos 2 6* 7 ), where Q y is the polar angle of the 
emitted photon with respect to the beam line. The invariant mass of the 
7T + 7r - pair is equal to the 77 mass: M vn =M v =547 .3 MeV. 

The KLOE detector consists of a large- volume drift chamber [3] (3.3 m length 
and 2 m radius), operated with a 90% helium-10% isobutane gas mixture, and 
a sampling calorimeter [4] made of lead and scintillating fibres. The calorimeter 
consists of a cylindrical barrel and two endcaps providing a solid angle coverage 
of 98%. A superconducting coil surrounds the entire detector and produces a 
solenoidal field 5=0.52 T. 

Tracks are reconstructed in the drift chamber with a momentum resolution of 
&(p±)/p±. < 0.4%. Clusters of energy deposits in the calorimeter are classified 
either as associated to charged tracks (charged pions, electrons or muons) or 
as isolated (photons, K^). Photon energies and arrival times are measured 

with resolutions of a E /E = 5.7%/ y/E(GeV) and a t = 54ps/ yj £(GeV) © 50 
ps; impact positions are measured with a resolution of a few centimetres. The 
readout granularity is 4.5 x 4.5 cm 2 in the plane transverse to the fibres, 
and is segmented in five layers along the particles direction. The trigger [5] 
is based on the detection of at least two energy deposits in the calorimeter 
above a threshold that ranges between 50 MeV in the barrel and 150 MeV in 
the endcaps. The higher machine background rates at small angle requires a 
higher threshold in the endcaps. 

Samples of simulated events are obtained using the GEANFI [6] code based 
on GEANT: event generators for any specific final state, including decay dy- 
namics and radiative corrections, are provided together with the detailed de- 
scription of the geometry and the response of each sub-detector. 

The decay chain — > 777 — > 7r + 7r~7 is searched for by selecting events with two 
tracks of opposite charge with a vertex near the e + e~ interaction point and 
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one prompt photon matching the missing energy and momentum obtained 
from the 7T + ir~ pair and the kinematic variables. The vertex is required 
to be inside a cylinder 20 cm long and 8 cm in radius, with axis parallel to 
the beam line, centred at the beam spot position. The polar angle 9 t of each 
track is required to satisfy 45° < 9 t <135°. A prompt photon is detected as 
an energy cluster not associated to any track, with time of flight T c i, distance 
from the beam spot position R c i, and energy E c i satisfying the condition \T c i — 
R c i/c\ < 5a t (E d ), where a t is the energy- dependent time resolution. In order 
to match the missing energy and momentum obtained from the ir + n~ pair 
with the photon kinematics, the angle ip between the direction of the missing 
momentum and the direction of the photon is required to be less than 0.15 
rad. 

Each track is extrapolated to the calorimeter and is required to be associated 
with a calorimeter cluster. A major source of background is due to radia- 
tive Bhabha events, e + e~— >e + e~ , ~f. Rejection of these events is based on the 
shower energy deposition in the calorimeter, on the time of flight, which is 
different for electrons and pions, and on kinematics. A likelihood estimator is 
constructed using the following information: the total energy of the cluster and 
the maximum energy release among the five planes of the calorimeter; the en- 
ergy release in the first and in the last fired calorimeter plane; and \T c i — L/c\, 
where L is the track length from the interaction point to the centroid of the 
cluster. The probability density function for each variable is obtained using 
samples of unambiguously identified pions from 7r + 7r~7r° and ii + n~ events. 
The separation between 7r + 7r~7 and e + e~7 events based on the values of the 
likelihood estimators L + and L„ for positive and negative particles is shown 
in Fig. 1. 

/i + /i~7 and residual e + e~7 events are rejected using the so called track-mass 
variable Mt- Mt is the particle mass computed assuming the <p decays to two 
particles of identical mass plus a photon. Mt is given by: 

\P4> ~ Pi ~n\ = E,/, - JpJ+ M~ T - yfp 2 2 + M T (1) 

where p\ and P2 are the three-momenta of the two pions and and p^ are the 
total centre-of-mass energy and momentum respectively. The track mass value 
is obtained from tracking information only and is very weakly correlated to 
the likelihood estimator value. The requirement 129 < Mt < 149 MeV selects 
7r + 7r~7 events; see Fig. 2. This cut, together with the cut on if) described above, 
ensures that background due to — > 7r + 7r~7r° events is negligible. 

The M n7T spectrum of the selected events ranges from 2771^=279 MeV to 
M0=1O19.5 MeV. Apart from the hypothetical signal, the physical processes 
which give 7r + 7r~7 final states are e + e~ — > ti + tx~ accompanied by ISR (initial 
state radiation) or FSR (final state radiation), — > f (980)7 with / (980) — > 
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Fig. 1. Scatter plot of the likelihood variable for positive (L+) and negative (L_) 
particles in arbitrary units. The line is the cut applied to select pion and muon pairs 
(above) from electron pairs (below). 
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Fig. 2. Distribution of the track mass variable My showing the separation between 
the peaks due to 7r + 7r~7 events (right peak) and to /i + /^~7 events (left peak). 

7r + 7r~ and (f> — > p ± 7r ± with p^ 1 — > 7r ± 7. Analyses of the full spectrum in differ- 
ent photon angular ranges are reported elsewhere [7,8] . Here we are interested 
in the M nn region around the t] mass where the signal is expected to be. The rj 
mass region (500 - 600 MeV) of the M vn spectrum is dominated by ISR events 
with the radiated photon mostly at small polar angle. To reduce the amount 
of such events while keeping a reasonable acceptance, we require the photon 
to be at large angle (45° < 9 1 < 135°). 

From Monte Carlo simulation we find the overall signal efficiency to be, 
e s =(16.6±0.2 stat ±0.4 S2/s t)%. The 2% systematic uncertainty is estimated by 




4 

3 
2 

1 



5 



comparing the data and Monte Carlo distributions of the variables Mt and 
ip. The overall rejection factors for the backgrounds range between order 10 4 
for /i + /i~7 and 10 6 for 7r + 7r~7r° and e + e~ r ). 

The expected M n7T distribution for a possible signal is a Gaussian with a 
resolution of 1.33 MeV. Analysis of the similar and abundant decay K s — > 
7r + 7r~ shows that the Monte Carlo correctly reproduces the observed mass 
distribution (see Ref. [6]). Figure 3 shows the measured M nn spectrum in the 
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Fig. 3. 7r + 7r~ mass spectrum between 500 and 600 MeV superimposed to the fitted 
background. The Gaussian is the expected signal shape in arbitrary units. 

region around the r] mass, in 1.2 MeV bins, together with the expected form 
for the signal. No evidence for the signal is observed. The curve superimposed 
to the data is the result of a fit to the M n7T spectrum over a much wider 
interval, from 410 MeV up to 1010 MeV with a function that describes all the 
appropriate physical processes: 



dN 
dNfZ 



'da(ISR) da(FSR + f ) da(pyr)' 



dikL 



x £ x e(M n7r ). (2) 



In Eq. 2 above, £ is the integrated luminosity of the sample analysed, e(M 7r7r ) 
is the selection efficiency as a function of M nn for the full range of the fit, and 
the da/d terms are the differential cross sections for the various processes 
which contribute to the background. The fit has 7 free paramaters and gives 
a x 2 value of 75 for 84 data points in the 500 - 600 MeV region and a value 
of 539 for 488 data points in the full mass range. 4 out of the 7 parameters 
describe the pion form factor according to Ref. [9]: M p , T p , a and /3; the other 
3 describe the scalar contribution according to Ref. [10]: M/ , q^kk and gf on7T . 
The result of this fit does not change if we remove the data points in the signal 
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region. We use this result as the estimate of the background magnitude in the 
following. 



In order to determine an upper limit for the branching ratio, we have repeated 
the fit by adding to the previously estimated background a signal component 
represented by a Gaussian with fixed mean and width of 547.3 and 1.33 MeV 
respectively, and free absolute normalisation. The fit returns a number of 
signal events N s = — 8 ±24, compatible with zero. The result does not depend 
on the choice of the fit interval and bin size. The probability distribution of N s 
has been checked by generating a large number of histograms according to the 
background distribution and fitting each of them to get N s . The results of this 
simple simulation show that N s , in the case of no signal, is indeed Gaussian 
distributed with a mean compatible with zero and a width of 24. The 90% 
confidence-level upper limit on the number of events is obtained using the 
tables in Ref. [11]. We find N s < 33. 

Alternatively, we have used a polynomial parametrisation of the background, 
obtained by fitting the sideband regions (500 - 540 MeV and 555 - 600 MeV) 
only. Applying the same procedure to get N s , we obtain N s = — 10 ± 24 and 
consequently N s < 31. The systematic uncertainty due to the parametrisation 
of the background is small and we use the largest value for the limit. 

The total number of 77's in the sample, N v , is evaluated counting the number of 
~^ VI events with 77 — > 3tt°. The efficiency for this channel is e(4> —>■ 777, t\ — ^ 
3n°) = 0.378±0.001 S f at ±0.008 S j /st , where the 2% systematic error is dominated 
by the uncertainty on the detection efficiency of low energy photons. Using 
the known branching fraction BR(?7 — > 37r°)=0.3251±0.0029 [12] we obtain: 

A(rw37r°) 7 
N„ = ~ '— — = 1.55 x 10 7 3 



with a systematic uncertainty due to the knowledge of the efficiencies and of 
the intermediate branching ratio of 2%. 

Taking the result N s < 33, the 90% confidence-level upper limit is 

BR{7] -> 7T+7T-) = — -. '— — - < 1.3 x 10~ 5 . (4) 

v ' ' N v e(r) -> ii+ix-) v ' 



This result is the best obtained to date and is 25 times more stringent than 
the previous best limit. 
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